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/define MIN_HI. 1 
/define MAX_HI 2001 
/define HI_INDEX 1 



/define MAX IDX 4000 



/* /define 


TISSDE */ 


/* 


/define C4 




/* 


/* #define 


RDX */ 


/* 


/* #define 


WG */ 


/* 


/* #define 


DYN */ 


/* 



Tissue-equivalent epoxy plastic 
C4 plastic explosive */ 
RDX sheet explosive */ 
Water Gel explosive */ 
40% dynamite stick */ 



/* new way of determining low */ 

/define zl .0247 

#define z2 .01492 

/define z3 .265 

/define z4 112.6 

/define z5 25.198 

/define z6 .6218 

/define z7 .265 

/* define substance parameters */ 

/if def WG 
/define cl 9.732 
/define c2 6.108 
/define c5 1.218 
/define K0 .547 
/define KL .961 
/endif 



/ifdef RDX 
/define cl 9.732 
/define c2 6.108 
/define c5 1.218 
/define K0 .65 
/define KL .86 
/endif 

/ifdef C4 
/define cl 9.732 
/define c2 6.108 
/define c5 1.218 
/define K0 .6522 
/define KL .87 
/endif 

/ifdef DYN 
/define cl 570.46 
/define c2 4.352 
/define c5 .304 
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\ /define K0 .522 
/define KL .765 
/endif 



/ifdef TISSUE 
/define cl 3798 
/define c2 3.8837 
/define c5 0.993 
/define KO .655 
/define KL .825 
/endif 

double bh (double km) ; 
double bh (double km) 
{ 

return (cl*pow( (km+c5) ,c2) ) ; 

} 

double Kref (double Hi, double Km, double kO) ; 
double Kref (double Hi, double Km, double kO) 
{ 



return ( ( (Hi+bh (Km) ) *ku*KL) / ( (bh (Km) *KL) + (Hi*k0) ) ) ; 



double alpha (double km) ; 
double alpha (double km) 

{ 

return ((zl+(z2*km)-(z2*z 3) )/(km*km)) ; 

} 

double beta (double km); 
double beta (double km) 
{ 

return( (z4+( (z6-km) * (z5/ (z6-z7) ) ) ) /km) ; 

> 

double newlow(double h, double km) ; 
double newlow (double h, double km) 
{ 

return (h* (1/ (km+ (alpha (km) * (h/ (h+beta (km) )))))); 

} 

double find_Km (double hi, double Kair, double kref) ; 
double find_Km (double hi, double Kair, double kref) 

{ .---"'* 

/* find the Km that approximates the desired Kref given high val,k0 */ 
int x,bitval; 

double lsbval,approx_kref ; 
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lsbval =0.8; 
bitval =04 

for (x=0;x<8;x++) 

{ bitval= (bitval«l) 1 1 ; 

lsbval = lsbval/ (double) 2 . 0 ; 

approxjcref - (Kref (hi , ( (double) . 1+ ( (double) bitval* lsbval) ) , Kair) ) 

if (approx_kref < kref) 

bitval=bitval& ( Oxf e) ; 



return ( ( (double) bitval*lsbval) +. 1) ; 

double findKm_Low (double hi, double low); 
double findKm_Low( double hi, double low) 

* /* find the Km that approximates the desired Low given high val,k0 */ 
int x, bitval; 
double lsbval , approx_low ; 

lsbval = 0.8; 
bitval = 0; 

for (x=0;x<8;x++) 

{ bitval- (bitval«i) 1 1; 

lsbval = lsbval/ (double) 2.0; 

approx_low = (Low(hi , ( (double) . 1+ ( (double) bitval* lsbval) ) ) ) ; 

if (approx_low < low) 

bitval=bitval&(0xfe) ; 

return ( ( (double) bitval* lsbval) + . 1) ; 



/* create the histogram */ 

for (hint = MIN_HI; hint < MAX_HI; hint += HI_INDEX) 

{ h = (double) hint; /* Get hi double value.*/ 

/* Set up the header values and the KIdx */ 
Hdr [ HI_VALUE ] = hint; 
KIdx = 0; 



\ 
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/* Get the hi and lo kref */ 

hi_kre.f = Kref(h, 0.29, kO) ; 

lo_kref = Kref(h, 0.8, kO) ; 
k=lo kref; 



lastl - -100.0; 
diffl - 1000.0; 
while (k<hi_kref) 

{ 

km=find_Km (h , kO , k) ; 
kr=Kref (h,km,kO) ; 
l=Low(h,km) ; 

if (((l-lastl)<diffl)&&(km>.29)) 
diffl = 1 - lastl; 

lastl = 1; 

if (h>800.0) 
{ 

k=k*1.04; 

k=1.01*k; /* 1% bins */ 

} 

/* do it again, but use diffl to find values */ 

k»lo kref; 

ta=fTnd_Km(h,k0,k) ; 

l=Low(h,km) ; 

findl=(int)l; 



/* adjust diffl to a power of 2 */ 
tdiffl=0; 

while ((1 « (tdiffl+1)) <= (int) diffl) 
tdiffl++; 

kxn=findKm_Low(h, (double) findl) ; 
k=Kref (h,km,k0) ; 

/* Save the minimum low and the scale factor */ 
Hdr [MIN_L0] = findl; 
Hdr [ LO_SCALE ] = tdiffl; 

while (k < hi_kref ) 
{ 

km=f indKm_Low(h, (double) findl) ; 
k=Kref (h,km,k0) ; 

/* Save the necessary information into the values */ 
KrefTab[KIdx] = (float) k; 



* 
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KIdx++; — v ' • 

/* increment low */ 
findl += (1 « tdif f 1) ; 

/* increment bin count */ 
bincnt+=l ; 

} 

/* Now we have the table, write out the header then the table */ 
Hdr [MAX_LO] = findl; 

bwritten = write (fhndl, (char *)Hdr, sizeof (int) *4) ; 
if (bwritten != (sizeof (int) * 4)) 

{ 

printf ("Error writing f ile\n") ; 
return (1) ; 

} 

/* Now write out the kref vector */ 

bwritten = write (fhndl, (char *) Kref Tab, sizeof (float) *KIdx) ; 

if (bwritten != (sizeof (float) *KIdx) ) 

{ 

printf ("Error writing file\n M ) ; 
return (1) ; 

} 

} 

/* output bin count */ 

printf ("Total Kref bin count :%ld\n",bincnt) ; 



Detection algorithm for above histogram 

I 
★ 

* Function: 

* DoBox 

* Descrition: 

* Process the box. 
* 

* Usage: 

* DoBox (x, y) 
* 

* Inputs: 

* x - int : the x coordinate of the candidate pixel 

* y - int : the y coordinate of the candidate pixel 
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* Outputs: 

* None * 

*/ 

static void DoBox (int x r int y) 
{ 

int tx, ty; 

double diffH, diffL, diffK; 
double kreflo,krefhi,krefavg; 
/* int tmp; */ 
double mindiff; 

Pixel *midpxl = &ScanLine[y] [x] ; 
Pixel *pxl; 

/* Average the values for this pixel */ 
AveragePixel (x, y) ; 

/* See if we need to do this pixel */ 
if (midpxl->avghia > 2000.0) 
r etur n ; 



/* 

* Calculate the min difference value (this is calculated by using 

* twice the expected noise as the difference value) . 

*/ 

mindiff = (10000*0/ (100. 0+midpxl->avghia) ) ; 

/* Now loop through the pixels doing the box */ 
for (ty = y - BORDER; ty <= (y + BORDER) ; ty++) 
{ 

/* Get the pixel */ 

pxl = &ScanLine[ty] [x - BORDER]; 

/* Loop through the x */ 

for (tx = x - BORDER; tx <= (x + BORDER) ; tx++, pxl++) 
{ 

/* See if we need to look at this pixel (edges are no-nos) */ 
if (pxl->sobel) 
continue ; 

/* Average this sucker */ 
AveragePixel (tx, ty) ; 

/* Now difference the Hi AIRS */ 
diffH = midpxl->avghia - pxl->avghia; 

/* Now threshold it */ 
if (diffH < mindiff) 
continue; 

/* Now difference the Lo AIRS */ 
diffL = midpxl->avgloa - pxl->avgloa; 



Page 
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/* Now threshold it */ ; ~ ^ 

if ((diffL < mindiff) || (diffL = 0*0)) 
continue ; 

kref lo=LookupKref (pxl->avghia,pxl->avgloa) ; 
krefhi=LookupKref (midpxl->avghia,midpxl->avgloa) ; 

dif fK = dif fH/dif fL; 

/* Key lookup algorithm 
* Histogram generation algorithm has been fit to this ratio 
*/ 

kref avg= (kref lo*. 8) + ( .2*krefhi) ; 

/* See if we need to histogram this point */ 
if ((diffK < (krefavg+(MinThreshold))) 

| | (diffK > (kref avg+(MaxThreshold) ) ) ) 

continue; 

midpxl->histval++ ; 

} 

} 

if (maxhit<midpxl->histval) 

maxhit=midpxl->histval ; 
if (midpxl->histval > fomThresh) 

f om += (midpxl->histval - fomThresh) ; 

if ( (midpxl->histval > 0) && (midpxl->histval <200)) 
histpix [ midpxl->histval ] ++ ; 



